Microscopic study of 4a-particle condensation with proper treatment of resonances 
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The 4q condensate state for ^^O is discussed with the THSR (Tohsaki-Horiuclii-Schuck-Ropke) 
wave function which has a-particle condensate character. Taking into account a proper treatment 
of resonances, it is found that the 4a THSR wave function yields a fourth 0^ state in the continuum 
above the 4Q-breakup threshold in addition to the three 0"*" states obtained in a previous analysis. 
It is shown that this fourth 0"'" ((O^)thsr) state has an analogous structure to the Hoyle state, 
since it has a very dilute density and a large component of q + ^^C(Oj) configuration. Furthermore, 
single-a motions are extracted from the microscopic 16-nucleon wave function, and the condensate 
fraction and momentum distribution of a particles are quantitatively discussed. It is found that for 
the (O4 )thsr state a large a-particle occupation probability concentrates on a single-a OS orbit and 
the a-particle momentum distribution has a 5-function-like peak at zero momentum, both indicating 
that the state has a strong 4a condensate character. It is argued that the (OJ)thsr state is the 
counterpart of the 0^ state which was obtained as the 4a condensate state in the previous 4a OCM 
(Orthogonality Condition Model) calculation, and therefore is likely to correspond to the Og" state 
observed at 15.1 MeV. 



I. INTRODUCTION 

The concept of a-clustering, is essential to understand 
the structure of light nuclei [HQ. In particular, ^^O is 
one of the nuclei which have the most extensively been 
studied with the cluster model as well as with the shell 
model. While the ground state is well understood to have 
a doubly-closed-shcll structure, the excited state, such as 
the second 0^ state, is known as one of "the mysterious 
0"*" states" whose excitation energy, 6.05 MeVjis too low 
to be explained in terms of the shell model [3[. It had 
also been mysterious that the second 0+ state forms a 
rotational band with low-lying 2+ and 4+ states. Be- 
sides those states, the third 0+ state at 12.05 MeV has 
been difficult to be understood from the shell model as- 
pect as well. One of the main aims of the cluster model 
approach has been to describe these states from a dif- 
ferent perspective and, indeed, the studies based on the 
a-cluster model have succeeded in reproducing many ex- 
perimental data. It was clarified that the and 0;^ 
states have aiS)+^^C{0+) and a{D)+^^C{2+) cluster 
structures, where the a particle rotates in an S and D 
wave around the ^^C core being in 0+ and 2+ states, re- 
spectively 0-01 • The great success also comes from hav- 
ing further revealed that almost all of the excited states 
up to ~ 14 MeV have considerable large a-i-^^C cluster 
components 0, [1]. 

Besides a-l-^^C clustering, 4a clustering can also be ex- 
pected to exist around the decay threshold into 4a par- 
ticles at 14.44 MeV, as suggested by Ikeda et al. Q- As 
an intriguing subject investigating the 4a cluster struc- 



ture, the presence of a 4a linear-chain state still remains 
an open question (sl-fioj. Another possibility are dilute 
gas-like states composed of weakly interacting 4a parti- 
cles They are analogues to the Hoyle state, i.e. the 
0^ state in ^^C, well described by a gas-like 3a-particle 
structure with dilute density [l2|. In good approxima- 
tion it forms a product state of the 3a particles in the 
lowest 05* orbit of the corresponding mean field poten- 
tial [ni[iM3- A striking fact is that the Hoyle state 
wave function obtained more than 30 years ago by solv- 
ing the 3a RGM (Resonating Group Method) equation 
of motion [l7| is to very good approximation identical to 
an antisymmctrizcd product type wave function [isj , the 
so-called THSR (Tohsaki, Horiuchi, Schuck and Ropke) 
wave function, introduced in Ref. jll| , and slightly mod- 
ification in Ref. [l^. Thus, the state has been proposed 
to be related to boson condensation of a particles in in- 
finite nuclear matter [l^-[2l| . The new interpretation of 
the Hoyle state as an a condensate has provoked many 
theoretical and experimental works on a-particle conden- 
sation phenomena in light nuclei |22h30| . It should be 
mentioned that also for A ^ An nuclei such as "'^^B and 
^^C several authors try to find states which are analogous 
to the Hoyle state fsillsst. 

The first investigation of the 4a-particlc condensate 



state was performed via the THSR ansatz in Ref. [11 



where two excited 0"*" states were obtained. In particular, 
the third 0+ ((Oj]')THSR) state has an excitation energy 
close to the 4a threshold and a considerably larger r.m.s. 
radius, i?ims — 3.9 fm, than the one of the ground state, 
i?rms ~ 2.7 fm, and therefore the (05]')thsr state was re- 
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garded as a candidate of the 4a condensate state. Based 
on this theoretical prediction, an experimental search 
for the 4a condensate state was undertaken via inelas- 
tic a scattering on ^^O. A new 0+ state was observed, as 
the fourth 0+ state, at 13.6 MeV with a relatively large 
a decay width of 0.6 MeV [H]. The angular distribu- 
tion of the cross section of the inelastic a scattering to 
the new state well agrees with the calculated one using 
the (0|^)thsr wave function. The relatively large de- 
cay width of 0.6 MeV is also consistently reproduced by 
the THSR wave function [35|. These results lead to the 
conclusion 113, HI] that the (Oj^)THSR state may be 
assigned to the new O"*" state at 13.6 MeV. Thus the 13.6 
MeV state was considered to be a candidate of the 4a 
condensate state [33 |. 

In order to confirm the above mentioned THSR de- 
scription of "'^^O, recently the present authors investigated 
the 4a condensate state via the 4a OCM (Orthogonal- 
ity Condition Model) [11] . The four-body (4a) problem 
with respect to the a-a relative motion was solved in a 
model space large enough to describe the dilute 4a gas- 
like configuration, as well as a+^^G cluster and shell- 
modcl-like ground state structures, using the Gaussian 
expansion method (37| . A one-to-one correspondence 
with the observed energy spectrum up to the 0^ state 
(see also FIG.|4]) was observed. It further was shown that 
the calculated 0^ ((Og )ocm) state has large occupation 
probability of 61 % for 4a's in the single OS" orbit and also 
a large overlap with the a-|-^^C*(IIoyle state) structure. 
These results indicate that the (O^)ocm state can be re- 
garded to be the 4a condensate state as the analogue to 
the Hoyle state, to be identified with the sixth experi- 
mental O''" state at 15.1 MeV. The well-known a-cluster 
structures of the Oj and 0^ states were simultaneously 
reproduced, in agreement with previous OCM calcula- 
tions 0, Q ■ It was also shown that the O4 and 0^ states 
both have a-|-^^C cluster structures, where the a cluster 
orbits in a higher nodal S'-wave and in a P-wave around 
the "'^^C core with = 0^ and lj~, respectively. These 
states had not been discussed so far in the literature. 

In this work, we confront our previous analysis with 
the THSR wave function for ^^O [11] with the above men- 
tioned results from the OCM approach. As we explained 
above, the 4a OCM results tell us that the Og state at 
15.1 MeV is a good candidate for the 4a condensate state. 
In Ref. [ll| using the THSR wave function only three 0+ 
states were obtained, including the ground state. This 
scarcity of 0+ states in the THSR approach is easily un- 
derstood by the fact that it allows only for two limit- 
ing configurations, that is a pure Slater determinant for 
B = b and a pure a-particle gas for S ^ 6 [H] . Asymp- 
totic configurations like ^^C(O^) -I- a are absent. They 
may be represented in a rough average way by the inter- 
mediate 0+ state. In any case the situation is such that 
it incited us to reanalyze the previous study of Ref. [ll| . 
The important question is whether or not the counterpart 
of the (O^)ocM state can also be found with the THSR 
ansatz, which should be a further state located above the 



(OJ)thsr state. In this respect, we need to recall that, in 
the previous THSR calculation [ll|, a proper resonance 
boundary condition was not imposed even for the states 
above the 4a threshold. Thus we could not obtain reli- 
able results for higher 0+ states than the (0;^)thsr state. 
The other excited state then only mocks up in an average 
way the states which have dominant a-|-"'^^C structure. 

The purpose of this paper is to report on results of 
a 4a THSR calculation with proper treatment of the 
resonances. We will see that this calculation gives us 
a new 0+ state as the (O4 )thsr state, slightly above the 
4a threshold. By comparing the wave functions of the 
(O4 )thsr state and the (O(i')oCM state in the 4a OCM 
calculation, it will be clarified that both states well cor- 
respond to each other. In particular, we will show that 
the (OJ)thsr state has a larger r.m.s. radius, a larger 
amount of condensate fraction, and a more sharpened 
peak around zero momentum for the a-particle momen- 
tum distribution, than the (0;^)thsr state. Those prop- 
erties are comparable to those of the (O^)ocm state. This 
indicates that the newly found (O4 )thsr state should be 
considered as the 4a condensate state and is likely to be 
assigned to the 0^ state observed at 15.1 MeV. 

The outline of the present paper is as follows: In Sec. 
Ill Ai we introduce the THSR wave function and the mi- 
croscopic Hamiltonian. The present approach to reso- 
nances is briefly explained in Sec. Ill Bl which is, together 
with the THSR ansatz, apphed to ^'^O in Sec HITll The 
obtained wave functions of ^®0 are analyzed in detail in 
Sects. IIITbI and IIII CI The quantities based on a-particle 
bosonic degree of freedom, such as the amount of conden- 
sation, are discussed in Sec. IIII D] Further discussions are 
presented in Sects. ilHEl and IIII Fl Sec.|lV]is devoted to 
the conclusion. 



II. FORMULATION 
A. THSR wave function 

In general, the na cluster model wave function can be 
presented in the following form: 



^nairi, ■ ■ ■ ,r4„) = A[xiRl,R2, ■ ■ ■ ,Rn)(t)ai4'ao ' ' >a„] 

(1) 

with A the antisymmetrizer and (jjai the intrinsic wave 
function of the i-th a-particle, which is taken as a Gaus- 
sian, 



Poii oc exp 



(2) 



i<fc<;<4 



The wave function x for the c.o.m. motion of the a's 
with Ri = j{ri^i + ri^2 + fi,3 + ri^^) is, of course, also 
chosen translationally invariant, that is it depends only 
on the relative coordinates Rij = Ri — Rj or on the 
corresponding Jacobi coordinates. The spin-isospin part 
in Eq. ([l} is not written out but supposed to be of scalar- 
isoscalar form. We will not mention it henceforth. 
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In 2001, Tohsaki, Horiuchi, Schuck, and Ropke 
(THSR) [11[ adopted, as the na condensate type wave 
function, the following ansatz for x in Eq- lU, 



THSR 
Ana 



{B; Ri, R2, • ■ • , Rn) 

MB;Ri-XG)---MB;Rn-XG), (3) 



with the total CO. m. coordinate = {Ri + - ■ ■ + Rn)/n 
and ipo{B; R - Xg) = cxp(-2(H - Xg)^/B^), that is a 
Gaussian with a large width parameter B which is of the 
nucleus' dimension. The product of n identical 05* wave 
functions reflects the boson condensate character. This 
feature is realized as long as the action of the antisym- 
metrizer in Eq. ([T]) is sufficiently weak. This type of con- 
densate wave function has been known, in the meantime, 
to have considerable success, notably with an accurate 
description of the Hoyle state, proposing it as the first 
of a series of excited states in na nuclei with a-particle 
condensate character. Those states can be considered as 
the precursors of macroscopic a-particlc condensation in 
nuclear matter at low density [19]. 

The internal part of the Hamiltonian we adopt in 
the present study of ^^O is composed of kinetic energy 
— ^V|, with nucleon mass M, the Coulomb force 
and the effective nuclear interaction: 



H = - 



^ 2M 

i=l 



VI -Tg 
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(2) 



E ^ 

i<j<k 



(3) 
ijk ' 

(4) 

where the c.o.m. kinetic energy Tg is subtracted. In 
this calculation, we adopted the same nucleon-nucleon 
force with two-body 
adopted in Ref. 
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(2) (3) 

and three-body V ■, terms as 
Both two-body and three-body 
forces consist of finite-range Gaussian functions with a 
parameter set named Fl given in Ref. [39| . 

The wave functions of quantum states in ^^O can then 
be expanded using the Aa THSR wave function, like 



[m] 



)$4a(i3^")) 



(5) 



where $4a(-B(")) is the 4a THSR wave function which 
has the form of 

(6) 

The discrete variational parameters iJ^™) represent the 
generator coordinate of the Hill- Wheeler ansatz. The 
expansion coefficients /fe(i?'''"'') and the corresponding 
eigenenergy Ek for the fc-th eigenstate are obtained by 
solving the following Hill- Wheeler equation. 



E 



<i>4a(S("'')) 

x/fe(s("')) = 0. (7) 



B. Extraction of the resonance wave function 

One of the aims of the present study is to discuss the 
states obtained by the THSR ansatz with proper treat- 
ment of resonances. The calculation of the resonance 
state in the bound state approximation is usually done 
by diagonalizing the Hamiltonian with the use of a fi- 
nite number of square-integrable basis wave functions. 
The calculated positive energy eigenstates which are lin- 
ear combinations of the basis wave functions are divided 
into resonance states and continuum states. 

There are many practical methods to carry out this 
division. One of the methods, found in Ref. [40] by some 
of the present authors, was shown to give consistent re- 
sults with the ones obtained with a well known approach 
to resonances such as the ACCC (Analytic Continuation 
of the Couphng Constant) method [l!], [i^- We apply 
the method to the present study, and following the liter- 
ature (40j , briefly explain it in this subsection. 

Let us first consider an attractive pseudo-potential V 
which is added to the original Hamiltonian iJ, yielding 



H'{5) = H + SxV, 



(8) 



where (5 is a coupling constant used to vary the strength of 
the pseudo-potential. We diagonalize this new Hamilto- 
nian H' [5) by using the same set of basis wave functions 
as used for H. As we increase the coupling constant i5 
from the physical value, (5 = 0, the eigenenergy of any 
resonance state decreases, and eventually the resonance 
state is transformed into a bound state. By contrast, con- 
tinuum states exhibit almost no change in their eigenen- 
ergies as 8 increases. 
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FIG. 1. Eigenenergies of the Hamiltonian ([8]) are schemati- 
cally shown as functions of S. There are three crossing points 
near which the mixing of continuum and resonance compo- 
nents occurs. (See the main text for a detailed explanation.) 
Figure adopted from Ref. [iol ]. 



This behavior of the eigenenergies as functions of 6 is 
shown schematically in FIG. [1] In this figure, the curve 
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labeled (b) around 6 = represents the eigenenergy of 
the resonance state denoted as i/'b (^) ■ It decreases as S in- 
creases. This curve (b) then crosses the curve (a) , which 
is fiat around 6 = and is the eigenenergy of a continuum 
state denoted as ipii{6). After the closest approach of the 
two curves (b) and (a), the curve labeled (b'), which is 
the smooth continuation of the curve (a), and hence is 
disconnected from the curve (b), represents the eigenen- 
ergy of the resonance state '4'h{6). This reminds us of the 
Landau-Zener level crossing scenario, see e.g. [11]. The 
eigenenergy of the continuum state ipai^)j which is first 
represented by the curve (a), continues with an almost 
flat shape, except near the crossing points. Eventually it 
is converted into the curve (a'). In the case of the res- 
onance state iphi6) at ^ = 0, we can consider that this 
state contains a negligible amount of admixture of con- 
tinuum state components, since it is located far from the 
crossing points with the continuum states. 

However, in the case of the resonance state ipd{6) cor- 
responding to the resonance state curve (d) in FIG. [U 
the state is at 5 = already involved in the crossing 
with the continuum state curve (c). Therefore, '>pd{6) at 
6 = should contain a sizable admixture of the contin- 
uum state. This is the usual situation for a resonance 
state with a broad width, and the resonance wave func- 
tion obtained by using the bound state approximation 
cannot be trusted, because of the sizable admixture of 
the continuum state. 

In order to overcome this difficulty, we note the util- 
ity of the continuum state wave functions obtained us- 
ing the bound state approximation in the neighborhood 
of the resonance state, since they should contain sizable 
admixtures of the resonance state component. In our 
present simple example considered in FIG. [TJ we notice 
the utility of the continuum state wave function ipci6). 
We can express both states, '0d(^) and V'c('5), as linear 
combinations of the pure resonance state 4'd(5) and the 
pure continuum state ((5): 

MS) = a^dW + /3*c(<5), MS) = -/^^dW + a^,{6). 

(9) 

The method for the extraction of the pure resonance com- 
ponent ^'d('5) which we proposed in Ref. [i^l is obtained 
by noticing that the extension of the density of a res- 
onance state is far smaller than the one of continuum 
states. More precisely, we can determine a and /3 from 
the requirement that the quantity 



{^di6)\J2{r,-XG)'\^di6)) 



be a minimum, where 



^d{6) = aMS)~ PAiS). 



(10) 



(11) 



Here A is the mass number. This requirement is satisfied 
by diagonalizing the operator ^ ^g)^ in the 

basis states ipd{6) and ipc{6), and by choosing (a, —f3) to 
be the eigenvector belonging to the smallest eigenvalue. 



When the resonance component is distributed among 
several energy eigenstates, we need to diagonalize the 
operator ~ ^g)^ in these energy eigenstates. 

When the energy eigenstates contain two resonance com- 
ponents, the diagonalization of the operator X]jLi(''i ~ 
Xg)"^ gives two small eigenvalues. Two resonance states 
are obtained by diagonalizing the Hamiltonian H'{6) 
with the two wave functions which give the two small- 
est eigenvalues of the operator ^f^iirj — XgY- 

We are now able to calculate approximate values of 
resonance energies along the resonance energy curve. If 
a point on the resonance state curve is not influenced by 
any crossing points with continuum state curves, we can 
regard the energy of the point to be just the resonance 
energy. On the other hand, if a point on the resonance 
state curve is influenced by some near-by crossing point, 
we extract the resonance wave function using the proce- 
dure explained above and then calculate the expectation 
value of the Hamiltonian H'{6) with the extracted wave 
function, which is just the resonance energy at this point. 
We apply this method for extracting the (OJ)thsr state 
with the use of the THSR wave function in Sec. IIII Al 



III. RESULTS AND DISCUSSIONS 

A. Extraction of the resonance components for the 
(OJ)thsr state 

In the Hamiltonian the explicit form of the long- 

(2) . 

range part of two-body nuclear force V^j is given by 



V 



-5[MeV] exp 



2.5[fm] 



with the value of Majorana parameter ' = 0.75 and 
Pjj the exchange operator of spatial coordinates. Sup- 

(2) 

pose that we here artificially vary the parameter M{ 
around the physical value 0.75, or equivalently, vary the 
coupling constant defined as 6 



(12) 



(2) 



0.75 around 



6 = 0. This corresponds to the situation explained in 
Sec. Ill B[ where a pseudo potential Vi^f^"^ is added to the 
original Hamiltonian, like 

H'{6)=H 



i<j 



(PS) 



(13) 



with 



(PS) 



5[MeV]exp[-(^^)'](l-P,,)- (14) 



The energy eigenvalues and eigenfunctions can then be 
obtained as (5-dependent solutions of the following Hill- 
Wheeler equation. 



^ U^B^M H{6) - EM) -fMB 



x/fc(i3(™'),<5) = 0. (15) 
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Here the fc-th eigenfunction is given as, 

^fcW = 5]/fe(i3'''"\'5)$4a(i3(™^). (16) 

rri 

In this calculation, for any values of 5., the fol- 
lowing 20 points of the generator coordinate B^™) = 

Y/62 + 2(i?^™^)2 are adopted: R^q"^ = 0.5 r^""^ with 
To = 1.196 and m = 1, • • ■ ,20. The size parameter of 
the a particle is fixed as 6 = 1.44 fm. 




FIG. 2. (color online). The energy eigenvalues of the Hill- 
Wheeler Eq. psp (solid curves) which are measured from the 
threshold (doted curve). The crosses denote the binding 
energies of the fourth 0"^ state (j>R{S) in Eq. (|17p . 




FIG. 3. (color online). The region of 0.0 < 5 < 0.05 of FIG. E] 
is enlarged. The adopted eigenstates in the superposition of 
Eq. (|17|l are denoted by open circles. The crosses denote the 
binding energies of the fourth 0''' state 4>r{5) in Eq. (|17p . 



In FIG.H the binding energies EkiS) - E'lJi, with 



the 4a threshold energy, are plotted up to 10 MeV in 
the range —0.01 < S < 0.16. The 4a threshold energy 
is calculated as El^ = 4£'„ -112 MeV, with £■„ the 
binding energy of the a particle, which does not depend 

(PS) 

on values of 6, since the pseudo potential V^j does 
not contribute to the binding energy of the a particle 
with the intrinsic wave function 0^ in Eq. ([2]). With the 
increase of S from left to right, we sec that three curves 
go down and pass through the 4a threshold one after 
another. The curves are considered to be trajectories 
corresponding to a resonance or bound state, which has a 
compact structure and is sensitive to the variation of the 
strength of the pseudo potential, as discussed in Sec. Ill Bl 
On the contrary, the other curves, which are relatively 
straight, i.e. insensitive to the variation, are considered 
to correspond to continuum states with unrestrictedly 
large spatial extension. 

We should here notice the important fact that the 
fourth 0+ state should be at around 3 MeV above thresh- 
old and that it has been missed in the previous calcula- 
tion of Ref. [ll|. With the decrease of S, the trajectory 
for the fourth 0+ state passes through the 4a threshold at 
S ~ 0.14, and then from S ~ 0.03 the trajectory quickly 
becomes invisible. This is because the fourth state has 
non-negligible decay width as a resonance. 

This is the very reason why the fourth 0"*" state could 
not be found in Ref. [ll|, while the first, second, and 
third 0+ states can without problem be assigned to the 
single eigenstates of the Hamiltonian, '^i{S = 0) = 5*1, 
^2{5 = 0) = ^2 and *3(J = 0) = ^-g, in Eq. dH]), 
respectively. 

However, the method explained in Sec. IIIBI enables 
us to extract the fourth state, which must exist in our 
approach around 3 MeV above the 4a threshold. The 
fourth 0"*" wave function can then be expressed as 

0^(5) = ^/?,.*^ (,5), (17) 

k 

where the coefhcients Dk are determined so that the cor- 
rect resonance wave function 4>r{5) should have a smaller 
radius than the continuum wave functions with spatially 
extended structures: (j>R{S) should satisfy the following 
relation against the variation of {D^}, 

, 16 

— (<^«(<5)| XcriMS)) = 0. (18) 

i=l 

In FIG. H where the region 0.0 < S < 0.05 in FIG. H 
is enlarged, the adopted eigenstates are marked by open 
circles at given 6 values, and the binding energy of the 
extracted resonance state is denoted by crosses. For 
example, at i5 = 0, we utilize in the superposition of 
Eq. pT|) the sixth, seventh and eighth states: ^6(<5 = 0), 
\i!j(S = 0) and <i/s{S = 0) obtained by solving the Hill- 
Wheeler equation (|15p . the binding energies of which are 
marked by three open circles, while the cross at S = 
denotes the binding energy of the extracted fourth 0^ 
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state. We can also see in this figure that binding ener- 
gies of the extracted resonance states denoted by crosses 
are reasonably on the trajectory for the fourth 0+ state 
in a wide range of S values, indicating the validity of this 
treatment of resonances. 

In all subsequent sections, as a matter of convenience, 
we refer to the wave function of the extracted fourth 
0+ state at (5 = as simply 4*4 instead of (f>R{S = 0). 
Thus, in the following, we represent the wave functions 
of (O^)thsr, (OJ)thsr, (0;[)thsr, and (0|)thsr states 
as ^1, \I/2, ^E's, and ^1/4, respectively. 



B. 



Features of the (0]^)thsr— (OJ)thsr states 



We present the calculated spectrum for the (O^)thsr- 
(O4 )thsr states in FIG. EJ together with the results of 
experiment and the 4a OCM calculation We can 
see that the OCM calculation gives a satisfactory one- 
to-one correspondence with the experimental spectrum. 
This can be understood by the fact that in the calculation 
the relative motions of a particles are solved in a huge 
model space, which is spanned by Gaussian basis func- 
tions [331. This allows us to represent the 4a-particle 
gas, a+^'^C clustering, and shell-model configurations, 
which are realized in the (Of}")ocM state, the (0^)ocm~ 
(05')ocM states, and the ground state (0]'")ocm, respec- 
tively. Thus, the (0,^")ocm state, which has the a con- 
densate character, is to be identified with the experimen- 
tal 0^ state at 15.1 MeV as the 4a condensate state. We 
should note that the observed widths of the experimental 
O4 , 0^ , and Og states also consistently correspond to the 
calculated ones (see Table |I| . In the same table, we can 
see that the M{EO) values for the (OJ)ocm, (O3 )ocm, 
and (O^)ocM states are in good agreement with the cor- 
responding experimental values. 

On the other hand, the present THSR ansatz gives 
only four 0+ states including the ground state. Let us 
now recall that in principle the THSR wave function of 
Eq. ([6]) allows for only two limiting configurations, a pure 
Slater determinant for B = b and a pure a-particle gas for 
B ^ b. Thus, as already mentioned, the scarcity of the 
0+ states comes from the fact that the THSR wave func- 
tions can represent the asymptotic configurations like the 
a-|-"'^^C clustering only in a rough average way. In par- 
ticular, the 0;^ and 0^ states which are considered to 
mainly include a{D)+^^C{2+) and a{P)+^^C{l^) clus- 
tering configurations, respectively, can hardly be repre- 
sented by this wave function, since it does not possess 
nonzero spin waves with respect to the a-a relative mo- 
tions. We should, however, note that the THSR wave 
functions are appropriate for describing states which 
dominantly have a gas- like structure of a particles weakly 
interacting in a relative 5-wave, as the a condensate. 
Also the ground state with essentially a Slater determi- 
nant configuration is within the possibilities of the THSR 
wave function. 

As mentioned above, we concluded in Ref. |39| that 
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4a threshold 



EXP. 4a OCM 4a THSR 

FIG. 4. Comparison of the 0"*" energy spectra between exper- 
iment, the 4q: OCM calculation [3y] and the present calcu- 
lation. Dotted line denotes the Aa threshold. Experimental 
data are taken from Ref [Zj ^nd from Ref. [s"" 
state. See Table U for the values. 



for the Ot 



the (0(i")ocM state, which might correspond to the exper- 
imental Og state at 15.1 MeV, may be identified as the 4a 
condensate state. It is therefore important to see whether 
or not the THSR ansatz can give a state corresponding 
to the (Og )ocM state and the 0^ state at 15.1 MeV state 
as well. On the other hand, as mentioned in Sec. U it 
was discussed in Ref. [131 that the (Oj^)THSR state could 
be assigned to the O4 state observed at 13.6 MeV, which 
corresponds to the (OJ)ocm state. This means that we 
should now examine the (O^)thsr state found in the pre- 
ceding subsection, as an important candidate for the 4a 
condensate state. 

In TableU the binding energies, r.m.s. radii, monopole 
matrix elements with the ground state M{EO), and 
a-decay widths of the (Of )thsr-(04 )thsr states are 
shown. The reason why the present numbers for the 
(0^)thsr-(OJ|')thsr states differ from those in Ref. [lH 
is due to the difference of the treatment of c.o.m. mo- 
tion in the THSR wave function. In the present work the 
c.o.m. motion is completely eliminated, while in Ref (ll| 
this is approximately done. The (0;j')THSR, state has a 
large r.m.s. radius of 4.2 fm and the (OJ)thsr state has 
even a larger one of 6.1 fm, which are comparable to the 
values for the (OJ)ocm and (O^)ocm states in the 4a 
OCM calculation, respectively. 

The (O4 )thsr state as well as the (O^)ocm state come 
quite a bit higher in energy than the observed 0^ state 
at 15.1 MeV. This probably entails that these states have 
the extremely large r.m.s. radii, 6.1 fm and 5.6 fm, re- 
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TABLE I. The binding energies E — E^^, r.m.s. radii Rrms, monopole matrix elements M(EO), and a decay widths F, in units 
of MeV, fm, fm^, and MeV, respectively, where £"4^ — AE^ denotes the 4q threshold energy, with E^ the binding energy of 
the a particle. The values of the previous 4a OCM calculation [sH] and of experiment [33,13] are also shown. 
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Rrms M{EO) 
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-14.37 
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- 0.05 
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0.185 


(04+ 


THSR 
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0.14 


(O5 )oCM 


2.08 
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1.0 r 


- 0.05 


0.66 
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0.166 



spectively. These two features possibly go together, since 
a state high up in the Coulomb barrier will necessarily 
have a quite more spatially extended structure. 

The M{EO) values of the (0+)thsr and (0|)thsr 
states well agree with those of the (04)0014 and (O^)ocm 
states, respectively. These indicate that the (O4 )thsr 
state corresponds to the (Og )ocm state, and hence to the 
15.1 MeV state. It will be discussed in Sec. lIII Fl to what 



extent 
state. 



the (O3 )thsr state approximates the (O4 ) 



4 JOCM 




15 



20 



FIG. 5. (color online). Nucleon density distributions of the 
(0^)thsr-(OJ)thsr states defined by Eq. (fT9)l . Note that 
they are multiplied by r^. 



The nucleon density distribution defined by, 

16 

gir) = {^k\-J2Si\r^-XG\-r)\^k), (19) 



16 . ^ 

for the (0^)thsr-(04 )thsr states are shown in FIG. [S] 
While the ground state has a compact saturation density, 
the (O^)thsr state has a very dilute density structure. 
This is a common feature of Hoyle-like states. Since the 
effect of the antisymmetrizer in the THSR wave function 



will be negligible in the state with dilute density [45[, 
this state is considered to have the product structure of 
loosely bound 4a particles, contained in the THSR wave 
function. 

The analogy to the Hoyle state can be discussed di- 
rectly in the ^^O system by calculating the overlap ampli- 
tude with the a+^^C(Oj) structure. The 4a condensate 
state must have a large overlap with this structure, since 
it is natural that if one a particle is knocked out from 
the 4a condensate, the remaining ^^C nucleus is in the 
Hoyle state with the 3a condensate structure. This can 
be seen in FIGS. |6l[a)-(d), where the following reduced 
width amplitudes (RWA's) for the (0+)thsr"(04 )thsr 
states are calculated. 



yL=oir) 



(20) 



with >I'(^^C) being in the ground or Hoyle state of ^^C ob- 
tained via THSR ansatz. For the ground state (0]^)thsr, 
the RWA in the channel a+^^C(0^) has a two- nodal 
behavior due to the Pauli principle. Accordingly, the 
(OJ)thsr- (0;]")thsr and (OJ)thsr states have RWA's 
with three, four and five nodes, respectively. The Hoyle 
state component is then the most largely included in the 
(OJ)thsr state (see FIG. E^d)), though the (O^)thsr 
state also has a relatively large Hoyle-state component 
(see FIG. WLc)). We show in FIG. [7] the RWA's for the 
(Oj)ocM state. We can see that the amplitudes in the 
ground state 0+^^0(0^) and the Hoyle state + ^^0(0^) 
channels, for the (Og )ocm and (O4 )thsr states, well cor- 
respond to each other. This indicates that the (O4 )thsr 
state is the counterpart of the (Og )ocm state and most 
appropriate for the 4a condensate state, rather than the 
(O3 )thsr state. On the contrary, the (O^)thsr and 
(OJ)thsr states dominantly have a large amplitude in 
the ground state channel a+^^C(0^), reflecting rather 
compact structures for these states. 

Based on the i?-matrix theory [i^, we can calculate 
the a-decay widths Fl for the (0;j')THSR and (O4 )thsr 
states, which are shown in Table HI with the following 
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FIG. 6. (color online). RWA's ryL^r) defined by Eq. (O for (a) the (0+)thsr, (b) (0+)THsa, (c) (0+)thsr, and (d) 
(O4 )thsr states in two channels a + ^^C(Oj) (dotted curve) and a + ^^C(Oj) (solid curve). 



fomulae: 



=2Pi(a)7i(a), 
ka 



F|(H + Gi(H' 

7i(a)=0i(a)7^(a), 



(21) 



7w(a) 



3r 
2^ 



where k, a and fi are the wave number of the relative 
motion, the channel radius, and the reduced mass, re- 
spectively, and Fl, Gl, and Phio) are the regular and 
irregular Coulomb wave functions and the corresponding 
penetration factor, respectively. The reduced width of 
el{a) is related to the RWA's as, 6l|(a) = ^yl{a). In 
this calculation, the decay energies are all taken as given 
by the experimental values: The excitation energy of the 
(O4 )thsr state is assumed to be the one of the observed 
0^ state at 15.1 MeV. Similarly the excitation energy of 
the (0|^)thsr state is assumed to be the one of the ob- 
served 64" state at 13.6 MeV. The total a-decay width 



of 0.14 MeV for the (OJ)thsr state is in good agree- 
ment with the corresponding experimental value of 0.17 
MeV. This total width only comes from the decay into 
the a -\- ^^C(O^) channel. The partial width decaying 
into the a + ^^C(O^) channel is completely suppressed, 
up to 10~^ eV, due to the very small value of the pen- 
etration factor PL=o{a) in Eqs. (PT|) . This is caused by 
the sufficiently small decay energy of 0.28 MeV into this 
channel, in spite of the large overlap between this state 
and a + ^^C(Oj) wave function with a certain channel 
radius, as seen in FIG. IH^d). In contrast, for the decay 
into the a-l- ^^C(O^) channel, large penetration is caused 
by the large decay energy of 7.9 MeV. Nevertheless, the 
small overlap of this state with the a + ^^C(O^) wave 
function suppresses the decay into this channel, and as 
a result, the relatively small width of 0.14 MeV is re- 
alized, indicating that this state is unusually longlived. 
The partial decay width into the a + ^^C{Qi) channel 
is calculated with a channel radius a = 12.0 fm. This 
mechanism also explains the small width of the (O^)ocM 
state, and hence the observed small width of the 0^ state 
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a+i2C(0i-) 
a+i-C(2^) 
a+i2C(4t) 
a+i2C(lT) 
a+i2C(3T) 
a+i2C(02) 




following squared overlap amplitudes between the wave 
functions ^fc(-B) and ^Pfe, 



20.0 



FIG. 7. (color online). RWA for the (Og )ocm calculated in 
Ref. m. 



at 15.1 MeV as well, which is discussed in Ref. t45[- In 
the present calculation, we neglect the decay into the 
other open channel a + ^^G{2f), since the THSR wave 
function of Eq. ^ practically cannot have nonzero spin 
components with respect to the a-a relative motions. 



C. Overlap with a single THSR wave function 



efc(i?o)= (^k'ifkiB)) , (A: = I,--- ,4). (25) 



Note that Ok{Ro) is written in terms of the parameter 
i?o, which is related to B like B = + 2B^. 




FIG. 8. (color online). Overlap amplitudes Qkijio) between 
the wave functions '^k and ^k{B) with k = I,-- - ,4, as a 
function of _Ro defined by Eq. (|25p. 



While the THSR wave function represents the physi- 
cally clear picture as the 4a condensation, the (O^)thsr, 
wave functions ^'fe with fc = 1, • • • ,4 of Eq. ^ arc ex- 
pressed as the superposed ones of the THSR wave func- 
tion of Eq. ^ with different B values. One might there- 
fore suspect that the superposed wave functions of the 
(O^)thsr states no longer keep the remarkable character 
of a single component wave function, retained in Eq. ^ . 
In order to clarify this point, we first construct orthonor- 
mal wave functions from the THSR wave functions as 
follows: 

*fc(B) =Cfe_i(B)Pfc„i$4a(5), (fc = l,---,4), (22) 
where Pk-i is projection operators defined by, 

Po = 1, (k = 1), 

fe-i 

Pfc_i = l-^|vl/0(vl/i|, (fc = 2,3,4), (23) 



1=1 



and Ck-i is a normalization constant. 



Ck-l{B) = (Pfe_i$4a(B) Pfc-l$4a(P 



-1/2 



(24) 



The wave function \E'fc(P) depends on the parameter 
B and is orthogonal to the wave functions ^'fc/ of the 
(0^/)thsr states with k' ^ fc. We then calculate the 



In FIG. [8] we show 6fe(i?o) of Eq. (fSSj) as a function of 
Po- Foi' any wave functions '^'^ (fc = 1, ■ • ■ ,4), Qk{Ro) 
(fc — 1, • ■ • ,4) takes values close to 100 % at optimum 
values of Rq. In particular, the maximum of 64 (Pq) 
amounts to 96 % at Pq = 6.5 fm. This means that the 
(O4 )thsr state is almost represented by a single THSR 
wave function parametrized by a large Pq value in the 
space orthogonal to the other (O^)thsr {k = 1,2,3) 
states. We should note that this feature is exactly veri- 
fied for the Hoyle state, strongly suppor ting the picture 
that the Hoyle state is a 3a condensate [Ij]. Therefore, 
the (04')thsr state shows the 4a condensate character, 
as does the THSR wave function $40 (Po) with a large 
Pq value. 

The (0;^)thsr state also gives a large squared overlap 
at a large Pq value, i.e. as a maximum 03(Po) = 0.98 
at Po = 4.0 fm. The r.m.s. radius of 4.2 fm of this state 
is comparable to the one of the Hoyle state, and there- 
fore this state also shows 4a condensate character. The 
(O^')thsr and (OJ)thsr states also have large squared 
overlaps, with maximal values of Oi(Po) = 0.98 and 
62 (Po) = 0.98 at Po = 1.2 fm and Po 2.5 fm, re- 
spectively. However, as shown in Tabled the r.m.s. radii 
of these states are rather small. In these states the an- 
tisymmetrizer strongly disturbs the product structure of 
4a particles, and therefore these rather compact states 
will not have a dilute gas-like structure of the 4a's. 
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D. Boson mapping of the THSR wave function 

The quantities such as occupation probabiUty of a 
single-a orbit and momentum distribution of a particles 
are very important information to discuss the nature of 
the a condensate. However, the wave functions ^'fc of the 
(Ofc )thsr states are not based on the degree of freedom 
of the a particles (boson) but on that of the nucleon 
(fermion), and therefore the various bosonic quantities 
cannot directly be given by the wave functions. In or- 
der to obtain them, it is necessary to extract the bosonic 
wave functions from the fermionic wave functions VPfe. 

The wave functions '•^k in Eq. ([S]) are rewritten as fol- 
lows: 



and 



Xk 



■^k=A{xm4>o.,---K,\, (26) 



(ov /fe(i?("))xir'''(s^"^ • • • , ^4) 



i=l 
3 



m i—1 ^ ^ 



(B(™))2 
2 i 



a I (27) 



where Xk is defined in terms of Jacobi coordinates: ~ 

R2 — Rl, ^2 = -^3 — {Rl + R2)/2, ^3 = i?4 — {Ri + 

R2 + Rs)/^. We here mention a few features of the wave 
hmction Eq. (|26|) . It can be rewritten in the form 



(28) 



where the following abbreviations arc used: (Prj 



and (, ^ 
defined by 



(4!)^ 
16!4! 



By using the norm kernel j47l |4 



(29) 



the orthonormal condition can be expressed as follows: 



(4!)' 



Skk' , 
(30) 

with d^l(d^l') EE d^j^d^^d^r^id^^d^'^dQ. The above re- 
lation suggests that, as discussed by several authors [T^ . 
[15I . I49I I , a normalized wave function in the bosonic space 
corresponding to the fermionic wave function '^k can be 
taken as 




dW/2(|,r)Xfe(^') 
J^'^'xk. (31) 



We should note that the above wave function is written 
in terms of the relative coordinates of the a particles, 
^, and docs not depend on the internal coordinates of 
a particles, which are integrated out in the norm kernel 

rnrri 

For simplicity, in Refs. p^, |49[ the following approxi- 
mate form has been taken as ipki^), 



^fxk 



VWxkWxk) 



(32) 



We also adopt this approximation in the following dis- 
cussions as one option. 

On the other hand, the square root of the norm kernel 
can be expanded in using the so-called exchange kernel 
K. MM, 



2 8 



(33) 



In dilute density states such as expressed by the THSR 
wave function, contributions from the exchange kernel K, 
can be expected to be small. We therefore approximate 
the above quantity, up to first order, as follows: First let 
us consider a projection operator A ^^jSl] onto the Pauli 
allowed space in the 4a system. Since the norm kernel 
TV is defined in the Pauli allowed space, the following re- 
lations are satisfied, M = 7VA = AJ\f — AM A, and hence 
J\f^^^ ~ AAf^^'^A. These properties should also be kept 
after the expansion of Eq. (p3| . which therefore should 
have the following form in the first order approximation: 



A/-1/2 ^ A (1 - /C)'/' A w A (^1 - i/C^ 



A 



(34) 



where for the second equality IC — 1 — JV is used. In 
diluted density states, to good approximation, the oper- 
ator A can only act on any pairs of a particles to exclude 
the Pauli-forbiddcn, OS', 01?, and IS relative states. We 
thus adapt A to the present calculation as follows. 



I 4 



(35) 



i<j s=QSfiD,lS 



where Us {Rij ) is the harmonic oscillator wave function in 
a OS, OD, or IS* state, denoted by s, for the a-a relative 
motion. Eventually, we can also adopt as another form 
the following approximate normalized wave function in 
the bosonic space, 



MO 



(A + A/')Xfe 



{A + N)xk {A+M)xk 



(36) 



We adopt the bosonic wave function -0^, in the forms 
of Eq. dnH) and Eq. for the present calculations of 

quantities such as the one-body density matrix in dis- 
cussing the amount of a condensation. Hereafter, we 
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refer to the ways of approximation for ipk of Eq. (|32)) and 
Eq. dSni) as APRl and APR2, respectively. 

By using the above a-particle bosonic wave function, 
the one-body density matrix for the a particle is defined 
following Refs. [H, [ll, 136,, i9. ^5l| , like 



Pair,r') = (?Afc|^ELift|V'fe), 



(37) 



p, = \SiR, -Xg- r')){S{R, - Xg - r) 
Then, by solving the following eigenvalue problem, 

dr'pa(r,r')(^(r') = (38) 



we can obtain the occupation probability /i and the corre- 
sponding single-a orbit Lp{r) as the eigenvalue and eigen- 
function, respectively. The occupation probability is 
labeled with the angular momentum L and the quantum 
number of a positive integer , like Ln^ : In this work, 
for a singie-a orbit with an angular momentum L, we de- 
note the largest occupation probability as Li {til = 1), 
the second largest as L2 {ul = 2), the third largest as 
^3 =3), etc. (see also Ref. [sa]). 
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FIG. 9. (color online). Radial parts of single-Q-particle orbits 
with L = Q and ul = 1 (5i) for the Hoyle state calculated 
with APRl (solid curve) and APR2 (dash-dotted curve) ap- 
proximations. See text for the definitions of APRl and APR2. 
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FIG. 10. (color online). Radial parts of single- a-particle or- 
bits with L = and ul ~ 1 (Si), for the (0]*')thsr (dotted 
curve) and (04')thsr (solid and dash-dotted curves) states. 
See text for the difference between the solid and dash-dotted 
curves. The values of jj, in the parentheses are the occupation 
probabilities into the corresponding orbits. 
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FIG. 11. (color online). Radial parts of single-a-particle 
orbits with L = and ul = 1 (Si), for the ground state 
(0]^)ocM (dotted curve) and the Hoyle state (OJ)ocm (solid 
curve) states, which are obtained from the 4a OCM calcula- 
tion [3a ]. The values of fj, in the parentheses are the occupa- 
tion probabilities into the corresponding orbits. 



In order to check if the above-mentioned two approxi- 
mations of APRl and APR2 work well, at least, in dilute 
density systems, we first apply the above formalism to the 
Hoyle state which has already been well known to have 
the dilute gas-like 3a structure. We show in FIG. |9]the 
single-a orbits for the Hoyle state, which are labeled as 
Si and obtained with the two ways, APRl (solid curve) 
and APR2 (dash-dotted curve). The microscopic wave 
function of the Hoyle state is given by solving the Hill- 
Wheeler equation with the THSR ansatz, where the same 
effective nuclear force is adopted as in the present study 



of ^^O. We can see that both curves are in good agree- 
ment with each other. This fact guarantees that the ap- 
proximations are realistic for the dilute density state. We 
see that the single-a 5* orbit has 05* nodal behavior with 
a very extended tail. This is the characteristic feature 
of the a condensate state [HI, [s^ , where the a particles 
occupy the lowest OS* orbit of a rather soft mean-field-like 
potential [l6j . 

As we discussed above for the ^^C case, also for -'^^O 
case it is important to calculate the occupation probabil- 
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FIG. 12. Occupation probabilities fi, which are labeled as L„j^ with L and n_t being the angular momentum and the quantum 
number of a positive integer, respectively, for (a) (0]'")thsr, (b) (OJ)thsr, (c) (0;^)thsr, and (d) (O4 )thsr states. The 
approximation of APRl is adopted to obtain the bosonic wave function. 



ity and the corresponding single-a-particle orbit, which 
are the crucial quantities for judging the amount of the 
a condensation. In FIG.[TOl the single-a S orbits for the 
(Oi')THSR (dotted curve) and (O4 )thsr (sohd and dash- 
dotted curves) states are shown. The solid and dash- 
dotted curves are obtained with the approximations of 
APRl and APR2, respectively, i.e. the ones of Ec^. ((32)) 
and Eq. ([55]) for the bosonic wave function of Eq. ([5T|). 
The dotted curve for the ground state is calculated with 
APRl . The occupation probabilities into the correspond- 
ing orbits are also shown in parentheses. 

The fact that both curves with APRl and APR2 for 
the (O4 )thsr well agree with each other again gives a 
support that both approximations are taken correctly. 
The approximation APR2 can be justified for the de- 
scription of dilute density states, so that the effect of the 



antisymmetrizer is weakened. 

We see in FIG. for the (O4 )thsr state the large 
single-a S orbit occupancy of 64 % and 68 % for the ap- 
proximations APRl and APR2, respectively, which has 
a nodeless behavior with a very extended tail. We can 
see m FIG. [m that the (Oe^)ocM state gives the similar 
value of 61 % of the occupation probability and the corre- 
sponding single-a S orbit also resembles the one for the 
(OJ)thsr state, concerning the OS" nodal behavior and 
the spatial extension. These results again give us sup- 
port that the (OJ)thsr state is to good approximation 
the counterpart of the (O^)ocm state described in our 
previous paper [s^. We should note that this behavior 
of the single-a orbit for this state is characteristic as an 
analogue to the Hoyle state (isj . 

On the other hand, for the ground state, the single-a 
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orbit has 25' nodal behavior with a strongly reduced tail, 
reflecting a compact density and hence the strong effect 
of the Pauli principle, for both approaches, THSR and 
OCM. 
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FIG. 13. (color online). Momentum distributions of the a 
particle defined by Eq. p9|l . (a) pa{k) and (b) k^pa{k), for 
the (0]^)thsr^(OJ)thsr states. The approximation of APRl 
is adopted to obtain the bosonic wave function. 
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FIG. 14. (color online). Momentum distributions of the nu- 
cleons defined by Eq. (|4H) . (a) g{k) and (b) k^g(k), for the 
(0^)thsr-(04 )thsr states. 



The distributions of the occupation probabilities into 
the single a orbits with APRl are shown for the 
(Of)THSR, (O^)thsr, (0;]")thsr, and (0|)thsr states 
in FIGS. [T^a)-fd). respectively. The orbits shown in 
FIG. [TU] correspond to the ones labeled as 5*1 for the 
(O^)thsr state (FIG. [HJa)) and the (0|)thsr state 
(FIG. [T^d)). We can thus see that the occupation prob- 
ability for the (O^)thsr state concentrates on the single 
OS orbit shown in FIG. [TU] with 64 %. This also means 
that the remaining 30 % to 40 % come from the occupa- 
tions into higher orbits, following the effect of the anti- 
symmetrization. On the other hand, the (0]^)thsr state 
has rather equal distributions of occupation probabilities 



labeled as 5*1, Di, and Gi. In particular, the largest con- 
tribution comes from the Gi orbit, reflecting more or less 
the doubly closed shell structure with the SU(3) config- 
uration of [4444] (A, ^) = (0,0). 

We should mention that the (0;^)thsr state and even 
the (02") THSR state also have considerable concentra- 
tion of occupation into a single-a S orbit as seen in 
FIGS.fTWb) and (c). The values of the occupation proba- 
bility are 40 % and 56 % for the (OJ)thsr and (O^)thsr 
states, respectively, to be compared with the value of 
64 % for the (O4 )thsr state. However, as mentioned in 
Sec, nil B1 we think that these states, just in a rough av- 
erage way, correspond to the observed Oj and Oj states 
with the a{S) + ^^C{Of) structures, respectively. It is 
more likely that the concentrations of the occupation 
into the single-a S orbits for these states disappear in 
a situation where the a-a relative motion is solved in a 
larger model space so as to include the a+^'^G cluster- 
ing configurations, as was done in the previous 4a OCM 
calculation. 

For the (0^)thsr-(04 )thsr states, the sums of the 
occupation probabilities for the ten orbits (^i. Pi, Di, 
Pi, Gi, S2, P2, D2, F2, and G2) are 0.93, 0.90, 0.90, and 
0.90, respectively. Due to antisymmetrization, there are 
about 10 percent missed from higher orbits other than 
those shown in FIG. 

The momentum distribution of a particle is also an 
important quantity to judge which state has a 4a con- 
densate nature. It can be obtained by calculating the 
following doubly Fourier transformation of the one-body 
density matrix pa{r,r') in Eq. ([ST]). 



Pa{k) 



1 



(27r)3 



(39) 



In FIG. [12] the above momentum distributions of the a 
particles are shown for the (Of )thsr~(0^)thsr states. 
As is consistent with the large occupation into the single 
a QS orbit, the (O4 )thsr state has the most prominent 
(5-function-like peak at around zero momentum, which is 
one of the typical characters of the a condensate state. 
According to the rather strong concentrations of the oc- 
cupation probability on the S orbits for the (O^)thsr and 
(OJ)thsr states, these states also have some enhance- 
ment around zero momentum. On the other hand, the 
(O^)thsr state docs not show such an enhancement. 

We show in FIG. [T3] the nuclcon momentum distribu- 
tions for the (of )thsr^(04 )thsr states, which can be 
defined, as in the a-particle case in Eqs. (|57]) and (|M]) . 
like 



kr' 



with 



e(r,r') = (vI/fe|^E-fift|*fe), 
g, = \5{r, -Xg- r')){5{r, - Xg - r)|. 



(40) 



(41) 
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The prominent peak of the a-particle momentum distri- 
bution seen at around zero momentum for the (OJ)thsr 
state disappears in the nucleon momentum distribution. 
This means that the nucleons, on average, have rather 
high momentum components even for the (04)thsr 
state. The nucleon momentum distribution corresponds 
to the Fourier transformation of the nucleon density dis- 
tribution, which is defined by Eq. (|19p (sec also FIG.[5|). 
Thus, this just reflects the diluteness of the nucleon den- 
sity and does not give the property of the condensate 
structure composed of the a particles. 



E. Size dependence of the occupation probabilities 

It is instructive to discuss the nuclear size dependence 
of the occupation probabilities of the single-a orbits in 
the 1^0(0+) state within the framework of the 4a THSR 
ansatz. The results are shown in FIG. [15] for the single- 
a orbits labeled as ^i, Pi, Di F^, and Gi. They are 
calculated based on the single 4a THSR wave function 
defined by Eq. which has the free parameter B trig- 
gering the size of the nucleus. Starting with the single 
4a THSR wave function, we extract the corresponding 
bosonic wave function with the approximation of APRl, 
and then calculate the one-body density matrix of the 
bosonic state, to obtain the occupation probabilities as 
the eigenvalues of the one-body density matrix. The ob- 
tained occupation probabilities then depend on the B- 
parametcr, i.e. the r.m.s. radius Rims of the nucleus. 




FIG. 15. (color online). Dependence of the occupation prob- 
abilities of the singie-a-particle orbits labeled as Si, Pi, Di , 
Fi, and Gi in the ^®O(0^) state on its r.m.s. radius -Rrms- The 
way of APRl is adopted to obtain the bosonic wave function. 



We can clearly see that with the increase of the r.m.s. 
radius i?rms the occupation probability comes to concen- 
trate on a single S orbit and the ones into the other orbits 
get rapidly suppressed. A similar behavior is also seen 



in the previous analysis for the '^^C(0"'") state based on 
the 3a OCM 15 1 and in infinite nuclear matter case j53| . 
where the condensate fraction is enhanced as the den- 
sity decreases. Here we notice that the r.m.s. radii 
of the (O^)thsr and (O4 )thsr states are calculated as 
-Rims = 2.5 fm and i?ims = 6.1 fm, respectively. The dis- 
tributions of the occupation probability at i?rms = 2.5 fm 
are close to the ones of the full solution of the (0]^)thsr 
state using the Hill- Wheeler equation, which are shown 
in FIG. [121 This is reasonable, since the single com- 
ponent THSR wave function, which coincides with the 
doubly closed shell wave function in the compact limit, 
as mentioned in Sec. Ill Al describes the (0]'')thsr state 
very well. We should recall that the squared overlap be- 
tween them is 0.98; as discussed in Sec. IIII CI However, 
the occupation probability of the Si orbit at i?rms = 6.1 
fm amounts to about 90 %, which is considerably larger 
than the value of 64 % for the (OJ)thsr state. The rea- 
son of this difference exists in the fact that the (O^)thsr 
state can no longer be represented to good approxima- 
tion by the form of the THSR wave function ^/ia{B) of 
Eq. ([BJ, but by the form of the wave function '^k=A{Ro)^ 
as precisely defined by Eq. (|^ . Let us now recall the dis- 
cussion made in Sec. IIII Cl where the latter wave function 
is constructed with the orthogonalization to the lower ex- 
cited states and the ground state. That is, we can think 
that the orthogonalization reduces the occupation prob- 
ability into a single S orbit to the value of 64 %, from 
about 90 %. In other words, we can say that the reduc- 
tion is caused by the effect of the antisymmetrization, 
which remains in the orthogonalization operator Pk-i=s, 
defined by Eq. 



F. The (OJ)thsr and (O^)thsr states in ^^O 

In the previous sections, we clarified that the (04^)thsr 
state apparently has a 4a condensate structure, and is 
likely to be the counterpart of the (Og )ocm state with 
the 4a condensate structure obtained in the 4a OCM cal- 
culation. We therefore consider that the (04')thsr state 
probably corresponds to the observed 0^ state at 15.1 
MeV, as does the (Og )ocm state. In this situation, we 
should keep in mind that the THSR wave function is not 
adapted to treat a-|-^^C configurations, since all a's are 
treated on equal footing. Therefore, the (O^)thsr and 
(0;^)thsr states obtained from the THSR wave function 
mock up in an incomplete way the a-|-^^C configurations 
in ^^O. It will be important in future work to include 
those configurations into the THSR description. 

On the other hand, the position of the (0;^)thsr state 
is rather close to the 4a threshold, and we see that the 
state also has a large amount of a condensate fraction, 
though the amount is not so large as that of the (O4 )thsr 
state. In this respect, the (0|^)thsr state could also be a 
candidate for the 4a condensate. However, as mentioned 
above, the result of the 4a OCM calculation indicates 
that a sizable fraction of the 4a condensation is only 



15 



included in the (O^)ocM state and not in the other five 0+ 
states. The explanation may again reside in the fact that 
we treat with THSR ansatz all a particles equally and 
then a relatively large r.m.s. radius value necessarily also 
leads to an a-gas-like state whereas in reality a large size 
also can be obtained from an a+^^C configuration. This 
points to the necessity to include a+^^C configurations 
into the THSR wave function. 

Let us mention that in Ref. [s^l the magnitude of the 
angular distribution of cross section in the (a, a') inelas- 
tic scattering to the observed state agrees very well 
with that calculated for the (0|^)thsr state, in spite of 
our present conclusion that the 0'^ state is not very simi- 
lar to the (O4 )thsr state. This agreement, nevertheless, 
seems rather natural, since as discussed in Refs. pvlls^l. 
the magnitude of the angular distribution in the inelas- 
tic scattering depends most sensitively on the size of the 
excited state and not so much on its internal structure. 



IV. CONCLUSION 

In this work, we reanalyzed the THSR wave function 
for the Aa system introduced previously by the last four 
authors of the present paper. 

Searching for states in the continuum we found an ad- 
ditional 0+ state. This (O^)thsr state was obtained at 
3 MeV above the 4a threshold and was shown to have a 
large r.m.s. radius of 6.1 fm. This finding could be done 
with the help of the proper treatment of resonances based 
on the ACCC method, which was developed by some of 
the present authors. We compared the wave function of 
the (O4 )thsr state with the one of the (O^)ocm state, 
which was previously obtained with the 4a OCM calcu- 
lation as the 4a condensate state. It was clarified that 
we could reasonably understand that the (OJ)thsr, state 
corresponds to the (O^)ocm state, and then is the most 
appropriate for the 4a condensate state, rather than the 
(O^)thsr state which had been considered to correspond 
to the condensate state in 11, 34, 35]. 

In the analyses of the THSR wave function, we showed 
that the (O4 )thsr state contains a large component of 
a+^^C{02) clustering in a way which is similar to the 
(O^I')oCM state. This indicates that both states corre- 
spond to each other and have a gas-like configuration of 
the 4a particles as the analogue to the Hoyle state, from 
the fact that the Hoyle state has a 3a-gas-like configura- 
tion as the 3a condensate state. The (O^)thsr state is 



then described by a single THSR wave function with a 
large value of the B parameter in a space orthogonal to 
the other three (0^)thsr-(03 )thsr states, though the 
(0]^ ) THSR" (O4) THSR states are in general expressed as 
the form of superposition of many THSR wave functions 
with respect to values of the B parameter. In order to 
further discuss quantitatively the amount of a conden- 
sate fraction, we extracted the boson degrees of freedom 
from the THSR wave function, which is microscopically 
described based on the nucleon degrees of freedom. It was 
done via somewhat approximate procedures which should 
be valid at low densities. We showed that the occupation 
probability for the (O4 )thsr state concentrates by a large 
amount on a single-a OS* orbit, like the (O^)ocm state. 
The shape of the orbits for both states is similar to each 
other and characteristic for Hoylc-like states. Further- 
more, the (OJ)thsr state is shown to have a strong peak 
around zero momentum in the a-particle momentum dis- 
tribution, indicating that the a particles are condensed 
around zero momentum. We also discussed the reason 
why the (O3 )thsr state, as well as the (OJ)thsr state, 
cannot be regarded as corresponding to the 4a conden- 
sate state, in spite of the fact that the states possess 
somewhat enhanced condensate-like features. 

AU these results strongly suggest that the (O4 )thsr 
state corresponds to the (O^)ocm state and is the most 
appropriately considered to be the 4a condensate state. 
This further gives us a strong support that the 4a con- 
densate state exists around the 4a breakup threshold and 
is very likely the observed 0^ state at 15.1 MeV. 

This conclusion will be supplemented with more ex- 
tended calculation than the present THSR description 
so as to incorporate a -I- ^^C clustering configurations. 
As concluded from comparison with OCM calculations, 
the a -|- ^^C components may then have a substantial 
contribution so that the condensate fractions will be re- 
duced. It is nevertheless expected that the (O^)thsr 
state remains an excellent approximation also in such an 
extended THSR description. 
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